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Charged Higgs bosons, appear in the particle spectrum of a general Two-Higgs Doublet 

Model (2HDM). We are concerned here with the case of Type II THDMs. In this context, the 
importance of singly produced charged Higgs bosons at future Linear Colliders (LCs) p| has 
been emphasised lately in several instances || g, |], g, [J. To begin with, it should be recalled 
that the detection of states would represent an unequivocal evidence of physics beyond the 
Standard Model (SM). Moreover, one may well face the following situation, as a result of the 
Large Hadron Collider (LHC) runs: only one light (below 130 GeV or so) neutral Higgs boson is 
found, h, and this is degenerate with the SM Higgs state. For example, this can happen over a 
large portion of the Minimal Supersymmetric Standard Model (MSSM) parameter space, in the so- 
called 'decoupling-limit', namely, when Ma 3> and for intermediate to high values of tan (3 (the 
two parameters that entirely define the Higgs sector of a 2HDM at tree- level). This also implies 
that the other MSSM Higgs states, H and H^, are similarly heavy (i.e., Ma ~ Mh ~ M H ±). 
Under these circumstances, one may have to wait for the advent of LCs, where precision tests 
of the Higgs sector can be performed, in order to fully clarify the nature of the Electroweak 
Symmetry Breaking (EWSB) mechanism. If the existence of a (light) Higgs state will have been 
proven at the LHC, then the most likely schedule at a LC will be to start running at a rather 
low energy (say, yfs = 350 or 500 GeV), where the corresponding Higgs production cross section 
(via e + e~ — » Z* — > Zh) is largest, as the latter proceeds via s-channel annihilation. At such 
energies, the heavier Higgs states may not be produced in the usual pair production channels 
j?]], either because below threshold (i.e., Ma + Mh,M h + + M H - > y/s) or since the intervening 
MSSM coupling in the decoupling limit becomes zero (e.g., in the ZZH vertex). Whereas in the 
neutral Higgs sector the heavy H and A resonances can always be accessed in 77 collisions (via 
77 — ► 'triangle loop' — ► H/A), this is not possible for the charged Higgs boson states, because of 
electromagnetic (EM) charge conservation. In the large tan (5 region, for neutral Higgs states, one 
could alternatively resort to the associate production mode e + e~ — ► bbH/A. The corresponding 
channel for a charged Higgs boson would be e + e~ — * btH + , which has an additional large mass 
in the final state (i.e., mt = 175 GeV). 

Hence, it becomes clear the importance of also studying production modes of charged Higgs 
bosons with only one such particles in the final state, in order to cover the mass range M H ±>y/s~/2 
at future LCs, where e + e~ — > H~ H + [|| falls short of a detectable cross section 0]. An analysis 
of the various single production modes was performed in Ref. 0, limitedly to their inclusive rates. 
(For similar studies of charged Higgs bosons produced in single modes in the case of ej and 77 
collisions, see Refs. @, 0.) There, it was shown that only two channels offer some chances of 
detection: 

e + e~ — ► t~v t H + , t + v t H~ (tree level), (1) 
e+e" -» W^H* (one loop). (2) 
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The former is relevant in the large tan (3 region, whereas the latter is important for the low one. 



As LEP2 data seem to prefer large values of tan /?, at least in the MSSM [11], we attempt here to 
devise a selection procedure that may help to extract process (Q) from the irreducible background^]. 
Since, as shown in Q, the production rates of process (|l|) are rather small in general over the 
mass region M H ± ^ yfsjl (for sake of illustration, we adopt here y/s = 500 and 1000 GeV), it 
is mandatory to resort to the main decay channel of heavy charged Higgs bosons, i.e., H + — > tb. 
Hence, the following processes are of relevance for the signal (S) and the irreducible background 

e + e~ — ► t~u t H + — ► T~v T tb (signal), (3) 
e + e~ — > it — > T~v T tb (background). (4) 

We require the emerging top to decay fully hadronically, i.e., t — > bW + — > 3 jets, whereas we have 
assumed t's to be tagged as narrow jets in their 'one-prong' hadronic decays: 

r± i^v T (12%), 

r± p ± (-> ttVK (26%), 

r ± af^vrVK (8%). 



Hence, the complete signature is: 



r-jet + ^ iss + 4j. (5) 



The simulation has been carried out at parton level. We have assumed a Type II 2HDM 
throughout with tan (3 = 40 and M H ± ranging between 160 and 660 GeV. For the signal, we have 
used the formulae of ]2| for the production process and the program described in [12] for the decay 
rates. For the backgrounds we have used the same code of Refs. [|b|, |l4|], also including non-it 
contributions, in which the two VF ± have been decayed appropriately^. All unstable particles 
entering the two processes faH* and W±) were finally generated off-shell (i.e., with the correct 
width). The integration over the final states has been performed numerically with the aid of 



VEGAS |T^] and Metropolis [19]. Finite calorimeter resolution has been emulated through a 
Gaussian smearing in transverse momentum, pr, with (ct(pt)/pt) 2 = (O^O/y 7 ^) 2 + (0.04) 2 for 
all jets. The corresponding missing transverse momentum, p™ lss , was reconstructed from the 
vector sum of the visible jet momenta after resolution smearing. A double tagging of 6-jets in the 



2 We defer a similar study of channel (|2|) to a forthcoming publication. 
3 Charged conjugated (c.c.) channels are assumed too throughout the paper. 



4 Also, we have verified that the background due to e + e — * bbZZ p{i}| , with one Z decaying hadronically and 
the other into two r's, one of which escaping detection, is negligible. Similarly, for the cases e+e" -> bbW ± H T and 
e + e~ — * bbH + H~ [[l6|, in which one boson decays to tau-nu and the other to light-quark pairs, and W~ +4j, with 
W~ — > r~v T , computed in [Q. 



3 



final state is implied. The non-running 6-quark mass adopted for both the kinematics and the 
Yukawa couplings was m& = 4.25 GeV. (The r-lepton was assumed to be massless throughout.) 
We neglect Initial State Radiation (ISR) and beamstrahlung effects, as we expect these to have a 
marginal impact on the relative behaviour of signal and background. 

We start our numerical investigation by comparing the LC rates for process (|l]) computed 
with the charged Higgs boson set on-shell (also refer to Ref. [g]) to those in which the latter is 
allowed to be off-shell. The corresponding curves are displayed in Fig. [|. No cuts are enforced 
here. At the 'threshold' point M H ± ~ y^s/2, one may notice that the two curves start departing. 
The effect is similar in size and shape at both energies considered. It is due to the finite width 
of the charged Higgs boson, which is of several GeV. By rewriting the Higgs propagator in the 
off-shell process as 



j + M H ± 



p 2 - M% ± + iM H ±r H ± V r t „, J <6) 

and taking the limit T^± — > 0, the Breit-Wigner in eq. @ becomes a representation of the 
Dirac delta function 5(p 2 — Mjj ± ) (apart from a factor ir) and the on-shell or Narrow Width 
Approximation (NWA) is recovered. For reference, in the same figure, we also show the rates 
obtained by using the two-body mode e + e~ — > H~H + followed by the decay H~ — > t~v t . 

The tails surviving the sharp drop at M H ± ~ are due to the diagrams that do not 

proceed via e + e~ — ► H~ H + — > t~v t H + and to the relative interference between the two sets of 
graphs. For reference, one should recall that the top-antitop background is at this stage (including 
the decay BRs yielding the four-body final state in (Q)) about 120 and 30 fb at = 500 and 
1000 GeV, respectively (at leading order). The S/B ratio is prohibitively large then, to start 
with, about 1:600(1:800) at M H ± sfsj% for = 500(1000) GeV. 

We now proceed our investigation by enforcing some selection cuts. Like in Ref. J2QJ , the 



Cambridge jet clustering algorithm |21j (see 22] for a comparative review of its properties) was 



enforced to isolate a five jet sample, here with y cut = 0.001, wherein the r-jet was treated on 
the same footing as the quark-jets. Similarly, both r- and quark-jets were required to pass the 
following cuts in energy and polar angle (hereafter, j represents a generic jet): 

Ej > 5 GeV, | cos%| < 0.995. (7) 

However, the former can be distinguished from the latter rather efficiently, thanks to very different 
sub-jet distributions (e.g., charged hadron multiplicities). Thus, one can apply a sequential 
and t mass reconstruction only to quark-jets, as follows: 

\Mjj - M w ±\ < 10 GeV, \M j:jj - m t \ < 15 GeV. (8) 

The cut in missing transverse momentum was: 

pf ss > 40 GeV. (9) 
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Finally, a very useful variable in distinguishing between signal and background is a transverse 
mass, Mt, constructed from the visible r-jet and the missing transverse momentum, i.e., 



M T = ^2p T T p™ ss {l - cos A(f>), (10) 

where Acj) is the relative azimuthal angle. In the case of the signal, the r-jets are heavily boosted 
relatively to the case of the background, as the charged Higgs masses considered here are much 
heavier than Myy±. By imposing 

M T > M w ± » 80 GeV, (11) 

the background is severely rejected while most of the signal survives. 

Now, recall that distributions of tt^ tracks coming from r's are sensitive to the polarisation 
state of the r-lepton and that in turn the spin/helicity states of r's coming from scalars and 



gauge bosons are different 24, In fact, the exploitation of this fact has already been 



proved to be very effective in the detection of H + — > t + v t signals in hadron-hadron collisions 



Basically, the key feature relevant to our purposes is the correlation between the polarisation 
state of the decaying boson and the energy sharing among the emerging pions. In fact, it is to 
be noted that the spin state of r's coming from H^- and W 1 * 1 -boson decays are opposite: i.e., 
H~ — ► TftVR and H + — » t^vl whereas W~ — > t£Dr and W + — > t^vl (neglecting leptonic mass 
effects, as we did here). Ultimately, this leads to a significantly harder momentum distribution of 
charged pions from r-decays for the .ff^-signal compared to the VF ± -background, which can then 
be exploited to increase S/B. This is true for the case of one-prong decays into both jr s and 
longitudinal vector mesons, while the transverse component of the latter dilutes the effect and 
must be somehow eliminated. This can be done inclusively, i.e., without having to identify the 
individual mesonic component of the one-prong hadronic topology. In doing so, we will closely 



follow Ref. 1 30] 



The mentioned transverse components of the signal as well as those of the background can 
adequately be suppressed by requiring that 80% of the r-jet (transverse) energy is carried away 
by the tt s, i.e.: 

R T = ^> 0.8. (12) 
Pt 

These arguments certainly applies in the region M H ± where the production cross section 

of process ([[]) is dominated by the two-body mode e + e~ — > H~H + followed by the decay H~ — > 



t~v t . Here, the enforcement of the constraint in (12) reduces significantly the background while 
costing very little to the signal. However, we are mostly interested in the complementary mass 
interval, M H ± ^ \ /r s/2, where pair production becomes irrelevant. Here, we have verified that the 
5 For the effects of r-polarisation in the neutral Higgs sector, see [E9|. 
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enforcement of the cut in (12) is harmless with respect to the effects on the S/B rates, as one can 
appreciate from Fig. ^, showing that the R T distributions for signal and background are basically 
identical. Therefore, we maintained the requirement (12) throughout the M H ± range considered 



here, for sake of consistency. 

A vigorous reduction of the background rates can however be obtained by enforcing all the 
other kinematic cuts described above: (|7|)-(|TT|). As intimated, Fig. [3| illustrates the strong impact 
of the constraint in transverse mass, by comparing the shape of the signal and background before 
the kinematic selection. The signal distributions are obtained at the points M H ± ~ y/s/2, at 
both centre-of-mass (CM) energies. 

The upper plots in Fig. || present the signal rates after the full kinematic selection has been 
enforced. (The background cross section is constant with M H ± as the above cuts do not depend 
on this parameter.) In the lower plots we display the significances (in black) of the signal rates, 
after 1 and 5 ab _1 of accumulated luminosity, C. It is clear that at this point neither evidence 
( ^ 3<7 ) nor discovery ( ^ 5<r ) of charged Higgs bosons is possible in the region M H ± > yfi/1, 
whereas for M H ± ^ y/s/2 the signals should be comfortably observed. 

At this point, one should recall that the charged Higgs bosons in the final state decay into 
visible objects, i.e., four quark-jets. Besides, recalling that the latter can efficiently be distin- 
guished from r-jets, one can look at the invariant mass of this multi-jet system. Fig. [5] shows this 
quantity. For the signal, it represents the reconstructed resonance of the charged Higgs boson, 
that we have originally generated at the point M H ± ~ y/s/2, for both CM energies. For the back- 
ground, it correspond to a non-resonant kinematic distribution. The width of the signal spectra 
is dominated by detector smearing effects and suggests that a further selection criterium can be 
enforced to enhance the S/B rates, e.g.: 

\M 4j - M H ±\ < 35 GeV. (13) 

The value of Mjj± entering eq. ( |l3| ) would be the central or fitted mass resonance of the region 
in M^j were an excess of the form seen in Fig. |5| will be established. 

The two red lines in Fig. |I]show the significances of the charged Higgs boson signals in presence 
of the constraint in (13), alongside those in (|7|)-(|i~2"D. By comparing these curves with the shape 



of the two-body H~H + cross section in Fig. |], one should expect to extend the reach in M H ± 
obtained from pair production of charged Higgs bosons and decays by about 20-30 GeV or so 
around and above the M H ± = \fs/2 point, at both CM energies, thanks to the contribution 
of single production. Typical signal rates at tan/3 = 40 in the threshold region would be 
6(30) events at y/s = 500 GeV and 2(10) at -/s = 1000 GeV, in correspondence of C = 1(5) 
ab _1 . Furthermore, recall that in the high tan/? region, where process (|l|) is of relevance, its 
production rates approximately scale like tan/3 2 , so that the higher this parameter the better the 
chances of isolating the signal discussed here. Finally, in our estimates so far, we have excluded 
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the efficiency of tagging the two 6-jets in the final state. According to Ref. |2C], the single 
6-tag efficiency is expected to be close to the value eb = 90%, so that our main conclusions should 
remain unchanged. 

We regard our findings as rather encouraging, especially considering the initial value of the 
S/B rates. Indeed, better selection procedures than those outlined here could be devised (see, 
e.g., |2(|), to improve further the discovery reach of heavy charged Higgs bosons, with masses 
M H ± > in a general 2HDM (including the MSSM). A realistic analysis exploiting more 



sophisticated simulations, based on the HERWIG Monte Carlo event generator [31, 32] interfaced 



to the TAUOLA package |33|] (see also [34|) for polarised r-decays, will be available in the near 
future 1531. 
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Figure 1: Total cross sections for process (|]) with the charged Higgs boson being on- and off-shell. 
No cuts have been enforced here. We also show the cross sections corresponding to e + e~ — ► H~H + 
production times the decay BR for H~ — > t~v t . 
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red). 13 
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Figure 5: Differential distribution in the invariant mass of the four quark-jets recoiling against 
the r-jet for the sum of processes (^) and (||) and for the former separately yielding the signature 
(|||), after the kinematic cuts in (]7|)-(|T^), including all decay BRs. 
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